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ABSTRACT 

At small-a:, the structure function Fl{x,Q‘^) is driven by the gluon content of the nucleon 
target and consequently it can unravel the underlying QGD dynamics in that region. In this 
work, one studies its behavior on photon virtuality at fixed energy within the color dipole 
formalism for models considering parton saturation effects. The reason is that they resum 
a wide class of higher-twist contributions, which have important influence in F^ description 
towards low Q^. It is shown that the geometric scaling property holds for the longitudinal cross 
section. Moreover, the effective anomalous dimension in scaling dipole cross sections can be 
investigated by studing both the turn over and the large regions of the recent experimental 
measurements. 


1 Introduction 


The longitudinal structure function Fl at high energies can be written in terms of the cross 
section for the absorption of longitudinally polarized photons. Its accurate measurements at 
low Bjorken variable x and/or photon virtuality would be helpful to constrain the underlying 
physics in that kinematical region. Namely, it is expected that it may discriminate between the 
leading twist predictions, which consider the collinear factorization and parton distributions 
determined from global fits, and the predictions from the saturation models which resum a 
class of higher twist contributions at small x. Along these lines, recently it has been claimed [1] 
that an accurate direct measurement could teach us the best way in which to use perturbative 
QGD for structure functions, since it gives an independent test of the gluon distribution at low 
X. It would be a much direct test of the success of different theories in QGD. The motivation 
for the following study is the recent extraction of Fl at fixed energy by HI Gollaboration at 
HERA [2]. 

Lets shortly summarize the present status of the perturbative QGD calculation using the 
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collinear factorization. The is well understood at high whereas very little is known about 
its behavior towards low and small x. Theoretically, we have that in the limit ^ 0 the 
structure function Fi has to vanish as reflecting the fact that the interaction of longitudinally 
polarized virtual photons has to vanish in the photoproduction limit. On the other hand, the 
leading twist DGLAP MRST [3] and CTEQ [4] global hts require the gluon distribution to 
be valence-like or negative at small x and low in order to describe the experimental data, 
leading to F^ being negative at the smallest x-Q^ bins. At that region, a comparison of the 
predictions at LO, NLO and NNLO using MRST partons have shown a poor description of 
the experimental results [5]. These results demonstrate the limitations of the applicability of 
perturbation theory and the necessity of resummation procedures. Nevertheless, in the global £t 
of the existing light-targets DIS data at LO, NLO and NNLO QCD approximations performed 
in Ref. [7], the high-twist contributions to the structure functions have been estimated. It was 
verified that these terms do not vanish up to NNLO and give important contributions at both 
small and large x regions. Therefore, a direct analysis of Fl could discriminate between leading 
twist and higher-twist resummations. 

There are also interesting and successful phenomenological analyses using the semihard 
factorization approach [8]. The longitudinal structure function is now obtained as a convolution 
of the structure function of the off-shell gluons with virtuality k± with the unintegrated gluon 
distribution. The data description is good, with the predictions being intermediate between 
the LO and NLO collinear calculations. The agreement is even improved if non-linear effects 
are included in the gluon distribution and/or NLO corrections to the approach are introduced. 
For the saturation effects, the GLLM unintegrated gluon distribution [9] has been considered, 
whereas an estimation for high order effects is taken into account by changing the hard scale 
in the strong coupling constant where = cQ^ (c 3> 1). The GLLM results for 

dipole cross section/unintegrated gluon function [9] are obtained by a global QGD analysis 
using the numerical solution for the Balitsky-Kovchegov nonlinear equation, which sums higher 
twists while preserving unitarity, and the short distance contribution is given by the DGLAP 
kernel. It has been shown that higher-twist contributions in the form of saturation effects result 
in better data description. This motivates an analysis on the dipole formalism [10], as k±- 
factorization approach and the color dipole approach are equivalent at the leading logarithmic 
approximation. 

An advantage in using color dipole models is that there are several analytical expressions for 
the dipole cross section, which describes the interaction of color dipoles with the (proton) target. 
This allows to obtain very fast calculations for the observables of interest and perform important 
qualitative studies on their behavior. In particular, it has been observed that the DESY- 
HERA data at small x and low can be successfully described with the help of saturation 
models. Moreover, the total cross section [11] and also the inclusive charm production [12] 
present the property of geometric scaling, which is a characteristic feature of the high density 
QGD approaches. The saturation (non-linear QGD) approaches are characterized by a typical 
scale, denoted the saturation scale Qgat(^) which is energy dependent and marks the 

transition between the linear (leading twist) perturbative QGD regime and saturation domain. 
The phenomenological models indicate that the saturation scale is smaller than 2 GeV^ at 
HERA and thus one expects that the signatures of the saturation effects become increasingly 
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evident in the region of small x and very low Q^. Furthermore, some of these approaches contain 
information of all orders in namely they resum higher twist contributions [13, 14], These 

corrections should be important at the low region, where the leading twist approaches would 
be in the limit of their applicability. These effects to Fl has been recently investigated within 
the dipole picture using saturation models at Ref. [15]. Early related analysis on longitudinal 
structure function within the dipole approach can be also found in Refs. [16, 17] and in Ref. 
[18] for a recent investigation using semihard approach. 

In this Letter, one analyzes the behavior of the longitudinal structure function Fl{x,Q‘^) 
at hxed energy as a function of Q^. One considers the color dipole formalism, focusing mainly 
on the saturation models. In particular, one investigates the analytical solutions for the color 
dipole cross sections presenting geometric scaling property. We also call attention that the 
longitudinal cross section at small-a; should exhibit scaling pattern as the inclusive case and 
this could be demonstrated from available experimental data. This Letter is organized as 
follows. In the next section, the theoretical description of Fl within the color dipole formalism 
is reviewed, introducing several implementations for the dipole cross section. In Sec. 3, one 
computes analytical expressions for Fl using scaling dipole cross section in the regions of large 

and in the turn-over region, which are useful in the determination of the effective anomalous 
dimension. In the last section, one presents the results for Fl at fixed energy and as a function of 
Q^, which can be contrasted with the recent experimental results. Moreover, one demonstrates 
that aL should naturally exhibit geometric scaling property and theoretical estimations are 
presented. Conclusions and a summary are also presented in last section. 


2 Structure function Fl and color dipole approach 


In the color dipole formalism, the deep inelastic scattering process can be seen as a succession 
in time of three factorisable subprocesses: i) the photon fluctuates in a quark-antiquark pair 
with transverse separation r ~ l/Q long after the interaction, ii) this color dipole interacts 
with the proton target, iii) the quark pair annihilates in a virtual photon. The longitudinal 
structure function is related to the longitudinal 7 *p cross section as Fl = —(j2 ^ ■ The 

latter is the overlap of the dipole cross section on the longitudinal photon wavefunction. The 
interaction is then factorized in the simple formulation [10], 

Fl{x,Q‘^) = J dz(fr\'^Liz, r; adipix, r), ( 1 ) 


where is the longitudinal momentum fraction of the quark in the color dipole, x = ^2 _^_q 2 

is equivalent to the Bjorken variable and provides an interpolation for the ^ 0 limit. The 
mass of the quark of flavour / is labeled as m/. The longitudinal photon wavefunctions 4/^ is 
determined from light cone perturbation theory and read as 




TT^ 


2:^(1 - zf K'^ {er) 


( 2 ) 
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where the auxiliary variable = z{l — z) + rrij and Kq is the modihed Bessel function. 
Finally, a summation over the quark flavors is understood in Eq. (1). It should be noticed 
that Fl goes to zero when ^ 0 at low x in the dipole picture since oc Q^. The dipole 
hadron cross section adip contains all information about the target and the strong interaction 
physics. There are several phenomenological implementations for this quantity and the main 
feature is to be able to match the soft (low and hard (large regimes in an unihed way. 

In the calculation presented here one considers analytical expressions for the dipole cross 
section, with particular interest for those ones presenting scaling behavior. Namely, one has 
<^dip oc {r'^QsatV for dipole sizes ^ 1/QLt where 7 is the effective anomalous dimension. 
In what follows one takes the phenomenological parameterizations: (a) Golec-Biernat-Wiisthoff 
model (GBW) [19], (b) Itakura-Iancu-Munier (IIM) model [20] and (c) Kharzeev-Kovchegov- 
Tuchin (KKT) [21]. In order to investigate the DGLAP evolution for small dipoles, the Bartels- 
Golec-Bienat-Kowalski (BGBK) model [22] is also considered. We call attention for the recent 
GLLM dipole cross section [9], which contains saturation corrections via solution of Balitsky- 
Kovchegov equation and DGLAP kernel for the linear regime. It is not included here as we 
are interested in analytical parameterizations. Let us summarize the main features of each 
parameterization. The GBW parameterization takes the eikonal-like form. 


CTdip {x, r^) 


o-Q 


1 — exp 





(3) 


where the parameters were obtained from a fit to the HERA data producing uo = 23.03 (29.12) 
mb, A = 0.288 (0.277) and xo = 3.04 ■ 10“^ (0.41 ■ 10“^) for a 3-flavor (4-flavor) analysis [19]. 
An additional parameter is the effective light quark mass, ruf = 0.14 GeV, which plays the role 
of a regulator at the photoproduction limit. The charm mass is set to be rric = 1.5 GeV. The 
GBW parameterization presents a geometric scaling form, adip oc f{r^Qlat)- For small dipoles 
< UQL it can be approximated by adip — <^o(’"^QLt/4), where the effective anomalous 
dimension is equal one, 7 = 1 . The BGBK model is the implementation of QGD evolution in 
the the dipole cross section, which depends on the gluon distribution as. 


^dip ^0 


1 — exp 


/ as{id^) xG {x^ lj?)\ 

1 , 3 ^;^ ) 


xG{x, Ql) = AgX (1 - xf-^ (4) 


where the initial condition is taken at Qo = 1 GeV^ and the hard scale is = G/r"^ + /ig. The 
phenomenological parameters are determined from a fit to small x HERA data. The function 
G{x, fi‘^) is evolved with the leading order DGLAP evolution equation for the gluon density. The 
improvement preserves the main features of the low-Q^ and transition regions, while providing 
QGD evolution in the large-Q^ domain. The effective anomalous dimension for the BGBK is 
expected to be close to the DGLAP anomalous dimension, 7dglap = 1- 


Although GBW and BGBK models are very successful in describing HERA data, their 
functional forms are only an approximation of the theoretical non-linear QGD approaches. An 
analytical expression for the dipole cross section can be obtained within the BFKL formalism 
and intense theoretical studies has been performed towards an understanding of the BFKL 
approach in the border of the saturation region [23, 24]. In particular, the dipole cross section 
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has been calculated in both LO and NLO BFKL approach in the geometric scaling region [25]. 
The IIM parameterization smoothly interpolates between the limiting behaviors analytically 
under control: the solution of the BFKL equation for small dipole sizes, r l/<5sat(a^), and 
the Levin-Tuchin law [26] for larger ones, r 3> 1/Qsa.t{x). A fit to the structure function 
F 2 {x,Q‘^) was performed in the kinematical range of interest. The IIM dipole cross section is 
parameterized as follows. 


^dip ( 2 ^) ^ 0^0 


'r^Q 


sat 


^ I In (2/rQgat) 
kXY 


0 (r - i?sat) + O-Q 


^ (brQaat) 


0 (i?sat - r) , (5) 


where i?sat {x) = 2 /( 5 sat and the expression for r( 5 sat(a^) > 2 (saturation region) has the correct 
functional form, as obtained either by solving the Balitsky-Kovchegov (BK) equation [27, 28], 
or from the theory of the Color Glass Condensate (CGC) [29]. The coefficients a and b are 
determined from the continuity conditions of the dipole cross section at rQsat = 2. The 
coefficients 7 sat = 0.63 (the BFKL anomalous dimension at the saturation border) and k = 9.9 
are fixed from their LO BFKL values. The parameters used here are ao = 27ri?p = 26 mb, 
A = 0.253, Xq = 0.267 x 10“"^ and Uq = 0.7. The light quark mass is set to be nif = 0.14 GeV 
in similar way as GBW model. The IIM parameterization presents scaling violation once the 
effective anomalous dimension depends also on the rapidity Y = ln(l/a;) for small size dipoles, 
7 {x, r) = 7 sat + For small dipoles having transverse size relatively close to i?sat the 

second term vanishes and the dipole cross section scales as adip oc (Tq 

The KKT parameterization has been proposed in order to describe hadron production in 
dAu collisions at forward and mid-rapidities. The phenomenological model for the quark and 
gluon dipole scattering amplitudes is inspired in the approximated analytical solutions of the 
BK equation for the saturation and color transparency regimes. The expression for the quark 
dipole-target forward scattering amplitude is given by [21]: 


^ X 'f(x,ry 

N Vsat 


Cdip {x, r) = (To - exp 
where the effective anomalous dimension 7 ( 0 ;, r) is parameterized as. 


7 (a^, r) = - 1 + 


K,{x, r) 


k{x, r) + J 2 k{x, r) -h 7C(3) 


K, {x, r) = 


ln|l/(r^Q;o)l 

(A/2)(y - y„) ■ 


( 6 ) 


( 7 ) 


where c is a free parameter and Q^^^iix) = A‘^x~^, with A = 0.6 GeV and A = 0.3. The initial 
saturation scale used in the expression for K,{x,r) is defined by = Qsat(Fo) with Yq being 
the lowest value of rapidity at which the low-x quantum evolution effects are essential. The 
form of the anomalous dimension is inspired by the analytical solutions to the BFKL equation. 
That is, in the limit r —a 0 at fixed x one recovers the anomalous dimension in the double 
logarithmic approximation 7 1 — In another limit of large Y with r fixed, Eq. 

(7) reduces to the expression of the anomalous dimension near the saddle point in the leading 
logarithmic approximation 7 i + uylp)- Therefore, 7 ( 0 ;, r) mimics the onset of the geometric 
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scaling region. In Ref. [21] the authors assume a characteristic value r 1/(2 /c_i_), where k± is 
the transverse momentum of the valence quark, and the anomalous dimension is approximated 
by 7 ( 0 ;, r) 7 ( 0 ;, 1/(4 A;^)). The parameters for KKT dipole cross section has been determined 
from the inclusive hadron production in deuteron-gold collisions at y/s = 200 GeV at RHIC 
(To = 0.6 and c = 4). However, it has been shown in Ref. [30] that those parameters are not 
able to describe the DESY-HERA data on structure functions. The main reason is that those 
parameters lead to a vary fast transition to saturation region and/or the saturation scale is 
somewhat larger (it goes to unity already at a; = 3 ■ 10“^) than GBW and IIM models. 

In order to obtain more reliable results with KKT model, in the analysis presented here 
the following modifications will take place. The saturation scale is taken from GBW model 
and the overall normalization is set as = {Nc/Cp) . With this procedure, the KKT 
model reproduces the GBW parameterization in the region where 7 ( 0 ;, r) = 1. Moreover, for 
computing k{x, r) one takes the identihcation = 1/Q^, Yq = 4.6 (which corresponds to 
X = 10“^) for the lowest value of rapidity at which quantum evolution is essential and c = 0.2. 
Putting together these ingredients, the KKT dipole cross section becomes similar to the BGBK 
dipole cross section. This fact seems to be consistent as BGBK has an effective anomalous 
dimension close to the anomalous dimension in the double logarithmic approximation for small 
size dipoles r — 0 . 

In next section, the analytic behavior of the longitudinal structure function at hxed energy 
as a function of is computed for dipole cross sections presenting scaling property. These 
results could help on a better understanding of the anomalous dimension. We will look at 
the large region, where small dipole configurations dominate, and also at the intermediate 
region of virtuality where saturation effects become important. 


3 Analytical results and effective anomalous dimension 


Let us start analyzing the large region. It is dominated by small size dipoles having trans¬ 
verse sizes lower than l/Qsat- In this region, the dipole cross section in the models presented 
above takes the approximate form adip oc (r^ QlatV■ Hence, at the limit r —0 GBW model 
gives 7 = 1, the IIM model gives 7 —1 (a rough pure scaling £t gives 7 ~ 0.84 for the range 
< 45 GeV^) and KKT provides 7 1 — Moreover, models based on BEKL evolu¬ 

tion equation present an anomalous dimension 7bfkl = 1/2. Using the scaling approximation 
for the dipole cross section it can be shown that. 


jdz d'^r 






4 a, 


N 

em4 ' c 


'Q 


l-7s 


sat 


TT 




E' 


(i-7s)r^(i + Y)r^(i-7s)r^(2-7s) 

47. r(2-27,)r(2 + 27,)(3-27,) 


( 8 ) 


where 7 = 1 — 7 ^ and the expression above is obtained in the massless limit. 

Inserting the result of expression Eq.( 8 ) into Eq. (1), the large behavior for the longi¬ 
tudinal structure function can be explicitly computed. For most of models 7 = 1, whereas for 
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BFKL approach 7 = 1/2. Thus, 7 ^ = 0 for GBW, BGBK and IIM (and probably for KKT) 
and 7 ^ = 1/2 for LO BFKL approaches. Therefore, one obtains. 




SiEeJQLW, for 7 = 1 (GBW), 

for 7 = 1/2 (BFKL) , 


( 9 ) 


The analytical results above give the theoretical expectation for the large region of the 
recent extraction of at hxed energy W^p = 276 GeV. As a; ~ Q'^/W^p, the saturation scale 
is running with virtuality and hence one has F^iW^pjQ^) cx ~ Q~^'^ for approaches 

having 7 = 1 at large virtualities. For the approaches based on LO BFKL evolution equation, 
one obtains Fi(lL/,p, oc This mild growth on virtuality is ruled out 

by data and the reason is that BFKL approach for structure functions should be valid just at 
intermediate virtualities. In fact, it is the transition for the double logarithmic evolution at 
small-a; and large Q^. For a pure scaling £t with 7 = 0.84 done in Ref. [], the behavior is still 
consistent and would give Fl{W^p,Q‘^) cx (q 2 ^i-o. 84 (i+a) ^ Therefore, the slope of the 

longitudinal structure function at fixed energy and large could help to theoretically single 
out the correct effective anomalous dimension. 


Now, let us look at the turn over region. This limit corresponds to the intermediate and 
low < 10 GeV^. This kinematical region is known to be in the geometric scaling region 
and the average dipole sizes are close to the saturation radius i?sat(a^)- The effective anomalous 
dimension is then given by the BFKL anomalous dimension at the saturation border 7 sat = 0.63. 
Therefore, at hxed x the longitudinal structure function behaves as Fl cx 
On the other hand, at fixed energy once again the saturation scale is running on and thus 
one has F7(hFyp,(5^) cx (q2^i-o.63(i+a) ^ Q‘’-2 . This corresponds to a mild growth on near 
to the turn over point. We can also try to estimate where the turn over point should be located. 
In order to do so, let us expand the dipole cross section around the value = 1. For 

the expansion, we take the functional form of GBW for sake of simplicity and consider just the 
two first term in its Taylor series, adip ~ (To( 1 — 2/e) + uo (r^ Qsat/^ e). These two terms are 
the main contribution to F^ in that region and the result can be obtained using Eq. ( 8 ). The 
next term would give the twist-4 contribution, which is a negative correction. For fixed cr, the 
longitudinal structure function will have the form. 


Fl{x, Q^) 


47r CKem 


QCem® f I Cjg 

-Uo 1 - - +- 

TT \ ej TT e 



( 10 ) 


where Therefore, the turn over region can be obtained by finding the point of 

maximum for the function above. It can be shown that it is located at cx Qsat(^)) with the 
proportionality constant being dependent on the prefactors of the dipole cross section expansion. 

Finally, it is timely to investigate the effective anomalous dimension as a function of Q^, tak¬ 
ing the association ^ 4 /for sufficiently large dipoles. It can be extracted from theoretical 
models by computing the following quantity, 7 efr = where the dipole scattering 

amplitude is related to the dipole cross section by J\f{x,r) = adip{x,r)/(7Q. The analytical 
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Figure 1 : (a) The effective anomalous dimension as a function of the quantity for 

different parameterizations. (b) The effective anomalous dimension as a function ofQ^. 
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results for GBW, IIM and KKT parameterizations read as, 

7T (rQs.t, V) 

(rQs^t, V) 


_ J 7sat 


ln(T-^QL,/4) 

kXY 


for rQsat > 2 , 

a \n{brQsa.t) exp [-a In^ {brQ sat)] : for rQsat < 2 . 


Cf 

Nr 


7' 


(r2 

< 5 slt)^exp 

1 


7 

47, 

1 — exp 

Cp 

Na 


7 

} 


7 = 7 (a:, Q^). 


( 11 ) 

( 12 ) 

( 13 ) 


In Fig. 1-a, the effective anomalous dimension is plotted as a function of It 

should be noticed that the saturation region is placed at r^Qsat/^ > 1- The GBW model gives 
a smooth transition from 7 —>■ 0 at the saturation region to 7 1 at dipole sizes r —0. The 

IIM model has a lower anomalous dimension than GBW in the saturation region, whereas they 
are similar at the scaling region rQsat ~ 2. The transition to 7 = 1 is slowest than GBW and a 
larger region is driven by 7 sat = 0.63 (1 < rQsat < 2). The discontinuity presented in the IIM 
plot is due to the step function for the interpolation between the different regions (saturation 
and color transparency). The behavior becomes steep in the limit r —> 0, strongly enhanced 
by the logarithmic dependence. The IIM parameterization violates scaling and for the plot we 
have considered x = 3 ■ 10“^ {Y = 8), which corresponds to saturation scaling being the unit. 
In fact, the effective anomalous dimension decreases with Y but at the same time is enhanced 
by the logarithmic contribution. The KKT parameterization presents also a slow transition to 
7 —>■ 1 , which depends on the virtuality through the identihcation r ^ l/Q"^ in the expression 
for 7 . In order to show this dependence, we have computed the effective anomalous dimension 
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for the typical values = 1 GeV^ and 100 GeV^ and fixed a; = 3 ■ 10 The obtained values 
are close to the saturation and BFKL anomalous dimension for a large span of dipole sizes. For 
Q2 = 1 GeV^, in the small dipole limit the KKT effective anomalous dimension goes to 0.68 
whereas for 100 GeV^ it goes to 0.84. In Fig. 1-a, for completeness we also show the DGLAP, 
saturation and BFKL values for the anomalous dimension. 

In Fig. 1-b is shown the effective anomalous dimension as a function of the virtuality 
Q^, using the average dipole size as r = 2/Q. The GBW parameterization presents a fast 
convergence for the DGLAP anomalous dimension at large Q^. On the other hand, both IIM 
and KKT parameterizations have a mild growth on virtuality and they converge for a value 
7 ~ 0.85 at large Q^. At intermediate virtualities, the GBW results have sizeable deviations 
in relation to IIM and KKT. At lower virtualities = 1 - 2 GeV^, the anomalous dimension 
reaches 0.63 for IIM and 1/2 for KKT. For completeness, it is shown also the KKT result using 
its original parameters. The effective anomalous dimension is very low in a large span of Q^, 
having an average value near the BFKL anomalous dimension. This fact explains the known 
impossibility of describing intermediate and large DESY-HERA data with the original KKT 
model. Its success in describing RHIG data on charged hadron spectra for d — Au is due to the 
low Pt 5 GeV range measured there and large mid-rapidity y ~ 0, where the saturation scale 
is not too large. It should be noticed that for the original KKT parameterization the saturation 
scale is extremely large, reaching about ^ 4 at x = 10“^. 

In short, the different anomalous dimension for each parameterization of the dipole cross 
section should account for a distinct behavior for the longitudinal structure function at hxed 
energy. Namely, the turn over point and the large region can provide valuable information. 
This is studied in details in next section. 


4 Results and Discussions 


Lets present the numerical results for the distinct parameterization at fixed energy W^p = 276 
GeV as a function of Q^, which are shown in Fig. 2. They are compared to the preliminary 
experimental data [2], namely the 99 min. bias data (circles) and 2000 shifted vertex data 
(squares). The total error for Fl is shown. The models give a reasonable data descriprion. The 
GBW parameterization gives the lower result and having the smaller turn over point, around 
= 20 GeV^. It presents also a less steep behavior at large virtualities. It has been shown 
it successfully describes the preliminary HI data in Ref. [2]. Therefore, a direct comparison 
among the several models is timely and necessary. The main shortcoming in the experimental 
measurements is large uncertainties leading to considerable systematic/statistic errors. More 
precise data and higher statistics would help to distinguish between different models. The IIM 
parameterization has a turn over point placed at at larger = 30 GeV^ and has a steeper fall 
at large Q^. Similar feature is present in KKT and BGBK models. This is due to the similar 
anomalous dimension for them in that region. The same particularity is present at the low 
Q^, where all parameterization have anomalous dimension near 7sat- The modihcation in KKT 
parameters presented in the previous section leads to numerical results having similarities with 
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Figure 2: The results for as a function of at fixed energy W^p = 276 GeV. The 

different numerical results correspond to distinct dipole cross section parameterizations (see 
text). 

to BGBK up to a different overall normalization. This can be contrasted with F 2 data and its 
universality could be tested. 

It is also timely to plot as a function of x at fixed values of as this is the usual 
way it is presented. It would be useful in order to figure out what are the differences in the 
x-behavior for different dipole models considered here. In Fig. 3, we show the results for the 
distinct dipole models in comparison with the experimental data [31]. The IIM and GBW 
models produce similar results for all Q^. Deviations are larger for BGBK at intermediate 
and very small x. The modified KKT model (with ad hoc parameters discussed in previous 
section) produces a consistent description, mostly at low . It can be improved by using a 
smaller overall normalization, which it would bring the theoretical curves close to experimental 
results. It is important compare these results from the color dipole approach against the more 
recent DGLAP QGD analysis. For doing so, we consider the last analysis on Fl using NLO and 
NNLO accuracy [32], which include the reduced cross section at high y and its earlier direct 
measurements. For comparison, we have chosen the IIM and BGBK dipole models. The reason 
is IIM being the more recent fit using saturation physics and it interpolates BFKL physics at 
large Q^. The BGBK model contains DGLAP evolution at large and saturation corrections 
at small-x and low Q^. For IIM, we include a multiplicative threshold factor (1 — xf to 
correctly describe the limit x —> 1. The results are shown in Fig. 4, for the sample virtualities 

= 2, 5, 20, 100 GeV^. There are strong deviations for low Q^, whereas the results are 
close at larger virtualities. At that region, the approaches are all convergent. The dipole 
approach gives numerical results which are between the LO and NLO/NNLO DGLAP analysis 
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Figure 3: The results for Fl{x, Q"^) as a function of x at fixed Q^. Data from HI Collaboration. 

and are consistent with the independent dipole model type £t in Ref. [33]. Concerning higher 
twist contributions, in Ref. [32] it is implemented via renormalon correction in the non-singlet 
sector. However, they basically remain the same behavior on x as the NLO/NNLO analysis. 
The deviation is larger at low-Q^, differing just by a distinct overall normalization. The dipole 
results are different from the renormalon correction as color dipoles models contain both higher- 
order corrections and higher twists. 

An interesting question concerning the longitudinal cross section is whether it exhibits 
geometric scaling property as the total inclusive cross section. Namely, the DIS cross section 
at high energies depends on x and exclusively via the scaling variable Tp = i-®- 

Q^) = (h)- Within the color dipole formalism and parton saturation approaches 

this property should be present at small-a: data. In order to address this issue, we consider two 
simple parameterizations for geometric scaling in the inclusive cross section and extrapolate the 
results for the longitudinal contribution. At low x and within the color dipole picture, the ratio 
between the longitudinal and total structure function can be roughly approximated by Fl/F 2 ~ 
2/11 [35]. For the functional form of the scaling function in DIS we take the results from 
Armesto-Wiedemann-Salgado in Ref. [34], where the high energy lepton-hadron, proton-nucleus 
and nucleus-nucleus collisions have been related through geometric scaling. The following 
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Figure 4: Comparison among color dipole results and recent NLO/NNLO DGLAP analysis. 

scaling curve for the photoabsortion cross section has been considered based on theoretical 
motivations: 


i^p) = (^0 [lE + F (0, /3) + In (/3)], P = a/Tp, (14) 

where je is the Euler constant and F (0, j3) the incomplete Gamma function. The parameters 
for the proton case were obtained from a £t to the small-x ep DESY-HERA data, producing 
a = 1.868, b = 0.746 and the overall normalization was fixed by do = 40.56 ph. The saturation 
scale is taken from the GBW parameterization. In Fig. 5 we plot (xp) as a function of 

Tp, where we have used the approximation a/, ~ (2/11) crtol- The functional behavior remains 
the same as for the total inclusive case. Experimentalists could study this property using the 
available data at small-a: for Fi and look for a geometric scaling pattern. 

The symmetric form of the scaling cross section has been also addressed in Ref. [36], where 
the symmetry between low and high region is related to the symmetry of the two-gluon- 
exchange dipole-dipole cross section. It has been assumed that in the dipole frame the energy 
evolution leads to the multiplication of partons and consequently to the appearance of the 
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Figure 5: The geometric scaling pattern for the longitudinal 7 *p cross section. The theoretical 
predictions are taken from scaling functions for the total inclusive cross section (see text). 


saturation scale When the virtuality is larger than the saturation scale, one 

has the usual DIS picture with the photon probing a set of independent partons. The growth 
rate of the parton densities is given by Isf, where A is such that the area of the proton 
is of order 1/A^ throughout the evolution. On the other hand, when the virtuality is close to 
saturation scale, the probability of multiple interactions between the photon and the proton 
becomes sizeable and a Glauber-like resummation is introduced and extrapolated by symmetry 
for virtualities lower than the saturation scale. Explicitly, the functional form for the scaling 
curve in the symmetric saturation model reads as. 





(l + log ./f^ 
i)Tp(l - logy 



}. 


for ^ 

for < Qsat • 


(15) 


where N (the normalization factor), a and b are free parameters hxed by DIS data. The result 
is plotted in Fig. 5 for the longitudinal case. It is similar to previous parameterization, with 
deviations for large Q^. This is due the limit of both functional forms at that limit. For 
completeness, we calculate the complete prediction from the GBW parameterization in order 
to check the approximation Fl/F 2 = 2/11 considered here. The approximation is reasonable 
in the intermediate region whereas it overestimates the ratio at low by about 30%. 


As a summary, we have analyzed the longitudinal structure function at hxed energy within 
the dipole formalism including parton saturation effects. These effects are important because 
they resume higher twist corrections to the process. In particular, we have shown the recent 
experimental data could help on investigating the effective anomalous dimension. The regions 
of interest are the large and the turn over point. We have computed expressions for the 
anomalous dimension for distinct analytical parameterizations for the dipole cross section and 
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show its role in the behavior. A precise measurement in the turn over point would allow 
distinguish among different models. It is shown that the longitudinal cross section should 
exhibit geometric scaling property as a natural extension as for the inclusive total cross section. 
This can be studied by experimentalists using the available measurements. Therefore, Fl is an 
outstanding observable testing both parton saturation and twist resummation and more precise 
data and/or more statistics is increasingly desirable. 
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